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The non-equilibrium dynamics of a system that is located in the vicinity of a quantum critical 
point is affected by the critical slowing down of order-parameter correlations with the potential for 
novel out-of-equilibrium universality. After a quantum quench, i.e. a sudden change of a parameter 
in the Hamiltonian such a system is expected to almost instantly fall out of equilibrium and undergo 
aging dynamics, i.e. dynamics that depends on the time passed since the quench. Investigating the 
quantum dynamics of a A^-component ip^-model coupled to an external bath, we determine this 
universal aging and demonstrate that the system undergoes a coarsening, governed by a critical 
exponent that is unrelated to the equilibrium exponents of the system. We analyze this behavior in 
the large-A limit, which is complementary to our earlier renormalization group analysis, allowing in 
particular the direct investigation of the order-parameter dynamics in the symmetry broken phase 
and at the upper critical dimension. By connecting the long time limit of fluctuations and response, 
we introduce a distribution function that shows that the system remains non-thermal and exhibits 
quantum coherence even on long timescales. 


I. INTRODUCTION 

Non-equilibrium behavior of interacting many-body 
systems is becoming a field of increasing importance in 
different areas of physics. It is largely driven by two main 
experimental advances: (i) the ability to bring a system 
into an out-of-equilibrium state in a controlled and re- 
producable manner and (ii) the ability to observe and 
follow the subsequent dynamics in real-time, i.e. on mi¬ 
croscopic timescales. Different tools have been developed 
for different experimental systems. Important examples 
are various cold-atom setups^^®, ultra-fast pump-probe 
spectroscopy of quantum materialsand heavy-ion 
collisions that explore the dynamics of the quark-gluon 
plasma^ These experiments clearly reveal that obser¬ 
vations far away from equilibrium can yield new insights. 

Non-equilibrium dynamics is often characterized by a 
lack of time-translational invariance. As a result, the 
fluctuation-dissipation theorem does typically not hold 
and memory effects occur, i.e. the dynamics depends on 
the initial state of the system. Major questions in non¬ 
equilibrium are (i) what are the effects of interactions; do 
non-linearities, aging and memory effects occur and how 
are they characterized^^^^^, (ii) what are the properties of 
transient metastable “pre-thermal” states^^^^^, and (hi) 
how does the system eventually reach a steady-state and 
is it given by a (generalized) thermal distribution^’^®“^°. 
Aging describes the phenomenon that correlation and re¬ 
sponse functions and depend on both 

time arguments t, t' and not just on their time differ¬ 
ence t — t' as is the equilibrium case, i.e. they depend 
the age of the system. These effects are well-known to 
occur in structural glasses, spin-glasses and disordered 
systems^^^^^. As these effects are related to spatial and 
temporal fluctuations it is often required to perform an 
analysis beyond the mean-field approximation. 

It is in general difficult to make quantitative predic¬ 


tions for the non-equilibrium dynamics of interacting 
quantum many-body systems beyond the mean-field ap¬ 
proximation^^- except in cases where exact solutions are 
available^®^^^ or when controlled numerical approaches 
are established such as for impurity models^®^®^ and for 
one-dimensional systems^^’®^^®^. Here, we exploit the 
presence of a quantum critical point to analytically solve 
for the universal dynamics of an interacting quantum 
model, which is an effective description of a number of 
experimental systems®^. We obtain solutions both for 
short as well as for long times. 

Universality close to (quantum) critical points is well 
established in equilibrium and follows from a divergence 
of the correlation length ^ and time U near critical- 
ity63,64 Observables such as the order parameter or cor¬ 
relation functions can be expressed in terms of univer¬ 
sal scaling functions with dimensionless arguments. As a 
function of, for example, the distance to the critical point 
they follow power-laws with universal critical exponents. 
In non-equilibrium, the correlation length is itself a func¬ 
tion of time. For example, if a parameter of a system is 
changed sufficiently slowly, the state and the correlation 
length of the system can adiabatically follow this change. 
Close to a critical point, however, adiabaticity demands 
that the parameter change must occur infinitesimally 
slowly due to the divergence of the correlation time. At 
any finite sweep rate, the system eventually falls out of 
equilibrium. In the Kibble-Zurek description of such a 
parameter sweep through a critical point, the correlation 
length is assumed to remain constant at this freeze-out 
length scale for the remainder of the sweep®^’®®. It then 
follows that the number of (topological) excitations de¬ 
pends on the rate via a universal scaling law that solely 
contains equilibrium critical exponents®^^^^ . Similarly, it 
follows that the long time approach to equilibrium of a 
system that is suddenly quenched close to a (quantum) 
critical point is governed by equilibrium exponents^^. We 
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confirm this behavior in our study as well. Kibble-Zurek 
scaling has been experimentally observed, in particular in 
the regime of slow sweeps^’^^^^®. In the opposite regime 
of fast parameter quenches also different scaling behavior 
has been reported^®. 

Recently, however, it was shown for the KZ regime of 
a gradual sweep, that the dynamics of a system does not 
completely freeze out once the system falls out of equi¬ 
librium^'. Instead, the correlation length ^(t) remains 
a function of time. If the system is located close to a 
critical point, where ^ diverges in equilibrium, the corre¬ 
lation length in fact diverges in a light-cone like fashion 
according to a power law 

^(t) oc (I) 

with a characteristic dynamic coarsening exponent ■ 
While the coarsening exponent is in general different from 
the dynamic critical exponent z that characterizes the 
quantum dynamics in equilibrium, it turns out that the 
two are equal in our approach. This scaling behavior (1) 
is also one of the results of our calculations for the case 
of a quantum quench. A crucial implication of it is that 
correlation and response functions obey a gener¬ 

alized scaling form, which incorporates the time depen¬ 
dence of ^(t). Using this scaling form, we are able to find 
an analytic solution of the universal dynamics of corre¬ 
lation and response functions, both at short and at long 
times. 

In this article we study the post-quench non¬ 
equilibrium dynamics of a dissipative Wcomponent 
model. Dissipation is introduced via coupling to an exter¬ 
nal bosonic bath with ohmic, sub-ohmic or super-ohmic 
spectrum. In equilibrium, the model exhibits a quantum 
critical point (QCP) that separates a disordered (symme¬ 
try unbroken) from an ordered (symmetry-broken) phase. 
The equilibrium properties of this model are discussed, 
for example, in Refs. 63, 78, and 79. We consider a sys¬ 
tem that is initially prepared in an equilibrium state away 
from the quantum critical point, where the correlation 
length ^ is finite. We discuss both the situation of an ini¬ 
tial state in the symmetry-broken and in the symmetry- 
unbroken phase. Then, at time t = 0 a parameter in 
the system is rapidly changed, i.e. quenched, to a value 
that corresponds to the quantum critical point in equilib¬ 
rium. Experimentally, this parameter can, for example, 
be pressure, strain, magnetic field or interaction strength. 

We are interested in the post-quench dynamics of the 
correlation and response functions as well as the time 
evolution of the order parameter itself. We obtain ana¬ 
lytical expressions for the universal part of the time de¬ 
pendence. Due to the rapid quench the correlation length 
first collapses to a non-universal value of the order of a 
microscopic length scale in the system. Since the system 
is located at the quantum critical point, however, it then 
recovers by critical coarsening in the light-cone fashion 
^(t) oc . An important implication of this power-law 
recovery of ^ is a universal prethermalized regime where 
the dynamics of the correlation function and the order 


parameter (j)(t) is described by the universal critical ex¬ 
ponent 

(j){t) oc . ( 2 ) 

For 9 > 0 follows that the order parameter increases 
after an initial collapse due to the fast spreading of cor¬ 
relations in the system, while for 9 < 0 the order de¬ 
creases. The duration of the prethermalization regime 
ix determined by the quench amplitude and increases for 
smaller quench amplitudes, i.e. if the initial state of the 
system was close to criticality. We note that while uni¬ 
versality is often associated with large length and time 
scales, it here arises from the critical light-cone recovery 
of the correlation length over time, similar to the coarsen¬ 
ing in heterogeneous systems®^. After the prethermalized 
regime at large times, we show that the system eventually 
relaxes to equilibrium quasi-adiabatically with a power- 
law that contains equilibrium critical exponents. Still, 
the prethermalization exponent 9 enters the long-time 
expressions as a universal pre-factor in the relaxation am¬ 
plitudes of the order parameter and the correlation func¬ 
tions. We note that the dominant slow dynamics is de¬ 
termined by the interactions that are present within the 
model and the coupling to the bath mainly ensures equi¬ 
libration in the long time limit. At long limes we can de¬ 
fine a time-dependent distribution function using 

the relation between response and correlation functions 
known as the fluctuation-dissipation theorem in equilib¬ 
rium. We Hnd that the deviation Sn from the equilibrium 
Bose-Einstein distribution is given by 

where 7 is the coupling to the bath, T the bath- 
temperature, denotes the equilibrium retarded 

Green’s funtion at the QCP and t is the time after the 
quench. As we show, this result implies that the dis¬ 
tribution function is non-thermal even in the long time 
regime: Sn decays only algebraically at large frequencies, 
while a thermal distribution decays exponentially. The 
approach to equilibrium over time is slow and described 
by a power-law oc . Finally, the fact that Sn changes 
sign and can become negative for z < 2 implies that the 
density matrix of the system is non-diagonal in the en¬ 
ergy basis and quantum coherence is present. 

The quench protocol to the critical point that we 
consider is inspired by the pioneering work on classical 
critical points by Janssen, Schaub and Schmittmann in 
Ref. 12 who analyzed the post-quench dynamics of a clas¬ 
sical theory in contact with an ohmic bath. The anal¬ 
ysis of the classical case was extended to colored noise 
in Ref. 81. Previously, we have reported results on the 
dissipative quantum (^^-model using a renormalization- 
group (RG) approach in non-equilibrium and calculated 
the new exponent 9 in an expansion ine = 4 — d — z 
close to the upper critical dimension®^. For a closed 
(/ 3 ^-model without an external bath, 9 was recently de¬ 
rived in Refs. 83 and 84 using a different non-equilibrium 
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renormalization-group method. In this paper, we com¬ 
plement our previous analysis by an expansion in 1/lV, 
where N is the number of components of the field. We 
also largely extend our previous study, which focused on 
a quench starting from the symmetry-unbroken phase, 
by a direct investigation of the order parameter dynam¬ 
ics (pit) starting from an initial state with non-zero 4>. 
We explicitly prove that the exponent 6 is independent 
of the initial state and the particular quench protocol. 
All previous investigations were limited to space dimen¬ 
sions below the upper critical dimension. We also discuss 
the behavior at the upper critical dimension due = 4 — z, 
where logarithmic corrections occur and 


(puclt) (X [\n{t/tj)Y"" (4) 


with tj being a microscopic time-scale that marks the 
beginning of the universal pre-thermalization regime. 

The remainder of the paper is organized as follows: in 
Sec. II, we introduce the model and the Hamiltonian, de¬ 
fine the quench protocol and the coupling to the external 
bath. In Sec. Ill, we present a non-equilibrium formu¬ 
lation of the large-A^ expansion and derive the large- 
saddle-point equations. To solve these equations self- 
consistently we employ a scaling form for the Green’s 
functions that we discuss in Sec. IV. We present the 
large-N solution for a quench to the quantum critical 
point starting from the disorder'd phase in Sec. V and 
starting from the ordered phase in Sec. VI. In Sec. V A, 
we relate the short-time scaling exponent 0 to the time- 
dependence of the self-energy, which can be captured by a 
time-dependent effective mass, and determine the short- 
time behavior of the Green’s functions. In Sec. VB, we 
elaborate on the long-time behavior of the Green’s func¬ 
tion, which enters into the central calculation in Sec. V G 
where we show the self-consistency of the solution, which 
fixes the value of the exponent 9. In Sec. VD we show 
that the system is characterized by a non-thermal dis¬ 
tribution function at long-times. In Secs. VE and VIG 
we discuss the dynamics at the upper critical dimension. 
In Sec. VI, we consider a quench that starts from the 
ordered phase, and find in Sec. VIA the scaling of the 
order parameter and the crossover time that separates 
the regimes of pre-thermal scaling dynamics from the 
adiabatic long-time dynamics that is characterized by 
equilibrium critical exponents as discussed in Sec. VIB. 
We finally summarize our main results and conclude in 
Sec. VII. In the appendices, we provide further details 
on the analytical calculations, in particular on the differ¬ 
ent short-time scales in the problem in Appendix A, on 
the derivation of the large-N equations in Appendix B, 
on the free post-quench Keldysh function in Appendix G 
and on the long-time limit of the full Green’s functions 
in Appendix D. 



Antiferromagnet Quantum Paramagnet 


Figure 1. (Color online) (a) Schematic description of the ip'^- 
model coupled to an external bath. The order parameter p 
experiences a potential landscape that depends on a set of 
parameters R. If the initial parameters Ri are such that the 
system is prepared in the symmetry broken state with finite 
order parameter (pi, we consider a sudden change of the pa¬ 
rameters from Ri to Rf that brings the system to the quantum 
critical point, where <(> = 0 in equilibrium, (b) One experi¬ 
mental realization of the (p'^ model for A = 3 is a quantum 
dimer model in contact to phonons. The quench can be per¬ 
formed by quickly rapidly the pressure. 

II. THE MODEL AND THE QUENCH 
PROTOCOL 

A. The ip^-model 

We consider the post-quench non-equilibrium dynam¬ 
ics of a <p"*^-model coupled to an external bath. The sys¬ 
tem is schematically depicted in Fig. 1(a). To be spe¬ 
cific, we also show one particular experimental realiza¬ 
tion, the dimer antiferromagnet realized, for example, in 
TlGuGls®® in Fig. 1(b). Other realizations of our theory 
are discussed below. The (^"^-model is described by the 
Hamiltonian 

(0 =) / (’T^ + (V¥>)^ + ro (() 

+ (5) 

Here, ip (x) = {(pi (x), ■ ■ ■, Pn (x)) is an V—component 
scalar field and tt (x) its canonically conjugated momen¬ 
tum with [pi (x) , TTp (a:')]_ = i6u'S {x — x'). In addition, 
we used the short hand notation f^ = f d‘^x. The non¬ 
equilibrium dynamics is a consequence of the explicit 
time dependence of either the bare mass rg (t), the in¬ 
teraction strength u (t) , or the externally applied field 
h (t), respectively. In the case of constant parameters 
the system is in equilibrium and undergoes at tempera¬ 
ture T — 0 and for h = 0 a quantum phase transition at 
a critical value rg ^ (u) of the bare mass rg. This transi¬ 
tion can be captured analytically by using an expansion 
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for small 1/7V. In fact, the coefficient 1/A^ in front of 
the term was introduced to allow for a well dehned 
N ^ oo limit. The large-fV approximation is comple¬ 
mentary to our earlier non-equilibrium renormalization 
group theory of Ref. 82. In particular, it allows for a 
rather straightforward analysis of the dynamics of a sys¬ 
tem that was initially in the symmetry broken phase. 


B. The quench protocol 

We consider the following protocols for the time de¬ 
pendence of 

R{t) = {roit),u{t) . ( 6 ) 

The system is initially prepared in the ground state with 
Ri = Wi, hi) somewhat away from the critical point 
located at Rc = (ro,c (u) ,u,0). Then we switch R(t) such 
that it approaches a final value R(t —>■ oo) = Rf that is 
closer to the critical point. The schematic phase dia¬ 
gram showing the parameter quenches that we consider 
in the following are shown in Fig. 2(a). The switching of 
R(t) takes place on a timescale r^. We consider the limit 
where Tg is of the order of the microscopic timescales of 
the system, i.e. the fastest time that is consistent with 
the validity of the model as an effective low energy 
description. This situation is depicted in Fig. 2(b). Keep¬ 
ing this in mind, we can safely perform the limit > 0 
in the calculation such that 

R(t)=R,+0(t)(Rf-R,), (7) 

with step-function 0 (t). 

After the quench the system instantly falls out of equi¬ 
librium. In the context of condensed matter systems, 
the generic situation is an open system that is coupled 
to external bath degrees of freedom. The bath, here con¬ 
veniently described by a set of harmonic oscillators, is 
assumed to stay at a constant temperature, T. In partic¬ 
ular, if T = 0 and Rf = R^ the system does indeed reach 
the quantum critical point in the limit t ^ oo. This is 
to be contrasted with the behavior of a closed system, 
where energy conservation alone implies that the system 
will not be in the ground state of the post-quench Hamil¬ 
tonian®^’®®. Even if the system thermalizes, a quench to¬ 
wards the critical point will heat up a closed system to a 
finite temperature. 


C. Coupling to an external bath 

If we include the coupling to the bath, the full Hamil¬ 
tonian reads 

( 8 ) 



Figure 2. (Color online) (a) Schematic equilibrium phase 
diagram of the ip^-model at T = 0 including the different 
quench paths that we consider in this article. We either 
approach the quantum critical point C from the symmetric 
side (path A —> C) or from the symmetry-broken side (path 
B ^ C). We also consider a quench in the presence of a 
magnetic field hi that induces a finite initial order parameter 
4>i. (b) Schematic plot of the quench protocol. We consider 
fast quenches where the switching time Ts from the initial pa¬ 
rameter set Ri to the final one Rf occurs on a microscopic 
timescale. 

Here, Hb describes the bath of harmonic oscillators and 
Hsb the linear coupling between system and bath: 

3 

Hsb^Y.^0 j ( 9 ) 

3 

In accordance with Eq. (7) the full Hamiltonian before 
and after the quench is then Hi = H {t < 0) and Hf — 
H {t > 0), respectively. 

In Eq. (9) flj are the frequencies of bath oscillators that 
couple with coupling constants cj. As the bath stays in 
equilibrium throughout, it is fully characterized in terms 
of the retarded function 

3 '3 

The imaginary part of ry (oj) is the spectral function of 
the bath. We consider a bath with spectral function 

Imry(a;) = (II) 

Here, the damping coefficient 7 determines the overall 
strength of the coupling to the bath and ujc is the ultra¬ 
violet cutoff of the bath spectrum. The exponent a char¬ 
acterizes the low-frequency behavior, a = 1 corresponds 
to an ohmic, a > 1 to super-ohmic and a < I to a 
sub-ohmic bath, respectively®^. The frequency depen¬ 
dence of Imry (w) for different a is depicted in Eig. 3(b). 


H{t)=H,{t)+Hb + Hsb. 
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A Kramers-Kronig transformation yields the real part of 
r] (w). We obtain 

6r] (w) = r? (w) - 77 ( 0 ) 

= 7 cot + isign (w)) |a;|“ ( 12 ) 

for frequencies small in magnitude compared to the cut 
off uJc- The zero-frequency value 77 (0) = depends 

explicitly on the value of the bath cut off. As we will 
see, 77 ( 0 ) merely shifts the location of the quantum crit¬ 
ical point of the system but does not affect the generic 
behavior near it. Below we will also need the analytic 
continuation 77 *^ (a;„) = 77 (0) -I- 6t]^ (w„) of 77 (w) to the 
Matsubara axis (w -I- iO'*' —)■ iuJn = i2mTT). It holds 

{Ujn) = , (13) 

sm 2 

again valid for frequencies that are small compared to 
the bath cut off: |u;„| uJc- 


D. Experimental realizations 

There are a number of experimental systems that 
can be effectively described by an Wcomponent ip"^- 
theory that we are considering. Realizations for 
A^ = 1 are the magnetic insulators CoNb 206 ®® and 
LiHoa;Yia;F 4 ®®, which can be described by the trans¬ 
verse field quantum Ising model. Dissipation in these 
systems arises via coupling to phonons. Realizations for 
N = 2 are systems near the superconductor-insulator 
(or superfluid-insulator) transition, Josephson junction 
arrays and quantum antiferromagnets in a magnetic 
field®^. Superfluid-insulator transitions in the Bose- 
Hubbard model have been experimentally investigated 
using cold-atom quantum gases confined to optical lat¬ 
tice potentials®®’®^ . Dissipation in those system has been 
engineered, for example, by coupling to other species®®. 
Dissipative nanowires near a transition to a supercon¬ 
ducting state®® and an ensemble of qubits in a photon 
cavity®"^’®® provide further realizations for N = 2. A 
realization for = 3 is the quantum dimer antiferro- 
magnet TlCuCla shown in Fig. 1(b), which can be driven 
across the quantum phase transition by either changing 
pressure or magnetic field®®. 


III. NON-EQUILIBRIUM FORMULATION OF 
THE LARGE-A EXPANSION 

The natural technique to treat the post-quench dynam¬ 
ics is the non-equilibrium many body formalism due to 
Schwinger®^, Kadanoff-Baym®®, and Keldysh®® (see also 
Ref. 100), where the field theory is placed on the three¬ 
time Keldysh contour shown in Fig. 3(a). A subtlety of 
our problem is that the initial state at t = 0 is itself 
a many-body state where (i) bath and system are en¬ 
tangled and (ii) interactions cannot be neglected. Thus, 


we cannot assume that we evolve from a non-interacting 
initial state and switch on the interactions adiabatically 
afterwards. How to modify the approach to the scenario 
where we prepare system and bath in an entangled inter¬ 
acting equilibrium state governed by Hi and where the 
subsequent time evolution is then determined hy Hf was 
discussed by Danielewicz®®® and Wagner®®®. In what fol¬ 
lows we will first summarize and then use this approach. 


A. Action within the three-time-contour formalism 


Let us consider an initial state of system and bath in 
equilibrium at temperature T. The state right before 
the quench at f = 0 is then characterized by the density 
matrix 


ft = 



(14) 


where /3 = 1/T and Z = is the pre-quench equi¬ 

librium partition function. In the limit T —>■ 0 it holds 
Pi = I'I'o) ('I'oIj where 1470 ) is the ground state of the cou¬ 
pled system and bath with Hamiltonian Hi. The subse¬ 
quent time evolution is governed hy Hf, i.e. the density 
matrix is given as 


p{t) = U{t,0)p,U{0,t), (15) 


where U {t, t') = exp (—7 {t — t') Hf) . The expectation 
value of an arbitrary operator (for specificity we consider 
our scalar field ip) is then 

{if (t)) =tT{p{t)ip) 

= ^ti:{e-^^'U{0,t)pU{t,0)). (16) 

For simplicity we suppress for the moment the component 
index I and spatial coordinate x of the field, i.e. ip (t) = 
Pi {x,t). The time evolution that enters the expectation 
value can be efficiently placed on the contour C shown in 
Fig. 3(a) with evolution operator such that 


{‘Pit)) = 




tr 


(Tee 


— i dsH (s) 


(17) 


The time argument of p denotes that the Schrodinger 
operator has to be inserted at time t. Te is the time 
ordering operator along the contour from i/3 ^ 0 ^ 
00 —>■ 0. In the denominator we used the fact that 
g-PHi _ jj jj g-PHi ^ analogy one can 

consider two-time correlation functions 


= -i {Tep {t) p {t')) 

= -*^tr (Tce-^^c (t) p (to) . ( 18 ) 

To determine these and other observables we consider the 
generating functional on the contour C = fc ^^) 


W[h] 
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r. 

ijS- 


(=) 

R{t) = Ri + 0{t){Rf — Ri) 


t 


lmri{Lo) / (b) 



Figure 3. (Color online) (a) Keldysh three-time contour, (b) 
Bath spectral function for different power law exponents cor¬ 
responding to the sub-ohmic case a < 1 , the ohmic case a = 1 
and the super-ohmic case a > 1 . 


with action 

5 [9?, X] = [cp] + Si, [X] + S,i, X]. (20) 

S'[<^,X] consists of the action of the system, the bath, 
and the coupling term. The individual terms are 

5's = ^ / {“ ’"0 

^ Jx,t&C ^ 

( 21 ) 

for the action of the system and 

Ssb = -Y.cJ X,(t).(p(t) (23) 

j Jx,t&e 

for the bath action and the system-bath coupling, respec¬ 
tively. 

We integrate out the bath variables at the expense of 
a bare propagator that is highly non-local in time: 

Go" {t, t') = - {dl + ro (t) - V^) d (t - t') 

-b A (t - t') • (24) 

The effects of the bath enter via the non-local self energy 
A{t — t') that is formally given as 

- c] {df + S{t- t'). (25) 

j 


The resulting action 

S[<f] = l [ ¥> (t) Go" (t, t') ¥> (!') 

^ Jx,t,t'ec 

- ^ J ^ . (26) 

depends only on the collective field (p and yields the gen¬ 
erating functional via 


W[h] = 




(27) 


In the usual Schwinger-Keldysh formalism it is conve¬ 
nient to place the four possible arrangements of the times 
t and t' on the forward and backward branch of the con¬ 
tour in a 2 X 2 matrix. In our case we have three branches 
of the contour which can be captured in terms of a 3 x 3 
matrix structure of the Green’s function^^^ 


G 


iG^ 

G< 

G< 

G> 

G^ 

G< 

G> 

G> 

G^ 


(28) 


The first matrix-element is up to a factor i the imaginary 
time Green’s function in equilibrium prior to the quench: 

(r -t') = - {TrPM (t) Pm (t')) , (29) 


where pm (t) = e'^^'pe~'^^P Tt is the time ordering 
operator along the vertical segment of the contour. The 
inverse of the bare Matsubara function is 

{gfy" (t t') = {dl - ro,i -b V^) 5 {t- t') 

-b (t - r), (30) 

where (r) is the Fourier transform of rj^ (w„) in 
Eq. (13). = G^(t — t') only depends on the 

difference between the two time variables. The Fourier 
transform of the bare Matsubara function yields 


(U^n) 


1 

-ujI - fo,i - -b St]^ (uJn)' 


(31) 


Here, we have introduced fo,i = ro^i — rj^(0) to account 
for a trivial shift of the bare mass due to the bath cou¬ 
pling. The functions 

(r, t') = i {pM (r) pn (t')) 

{T,t') = i{pH {t')pM (t)) (32) 

describe correlations across the quench, where pn (t) = 
is given in the Heisenberg representation 
after the quench. The remaining 2x2 block with 

G"^ = -i {TtpH {t) PH {t ')), 

G'^ = -i (ftPH it) PH (i')) . 

G> = -i {pH {t') PH it)) , 

G< = -i {pH {t) PH (<')) 


( 33 ) 
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makes up the usual Keldysh matrix where T* and Tt re¬ 
fer to the time ordering and anti-time ordering operator 
along the horizontal real-time branch of the contour C, re¬ 
spectively. Of those 9 Green’s functions that occur in the 
3x3 matrix structure, only a few contain truly indepen¬ 
dent information. By transforming Q with the rotation 
matrix 


L = 


-1 0 
0 ^ 



(34) 


gf (t, t') = {t — t') (both before and after the quench) 
and can be Fourier transformed to yield 


3/(w) 


1 

- foj - q'^ + Srj(uj) ’ 


(40) 


where foj = tqj — Notice, that this form is only 

correct for the bare retarded function, i.e. for u = 0. 
As we will see, many-body interactions couple response 
and correlation functions and lead to aging effects in the 
retarded function as well. 


one can benefit from this redundancy among the Green’s 

functions and obtain Large-N equations 


/ G^ 

-V2G< 

0 \ 


-V2G> 

G^ 

G« . 

(35) 

V 0 

GA 

0 / 



The scalar field forms a three component vector 
(<PM 7 , <p|i')with so called classical and quantum com¬ 

ponents 

(t) = (v’h (t) ± Vh (t)) > ( 36 ) 

in addition to the Matsubara component. The two rel¬ 
evant Green’s functions on the forward and backward 
branch of the contour are the retarded Green’s function 

= -iO {t - t') {Yph {t) (37) 

and the Keldysh function 

{t, t') = -i {[ifH (t) , (fH (i')]+) . (38) 

where [A, B]± = AB ± BA. The retarded function mea¬ 
sures the response of the order parameter at time t caused 
by an external field that couples to it and that acted at 
time t' < t after the quench, while the advanced function 

G^ = lO it' - t) {[g^H it) , (39) 

contains the same information as the retarded function: 
G"^ = G^{t',t). In distinction, the Keldysh func¬ 

tion is a measure of the strength of correlations between 
order parameter configurations at distinct times after the 
quench. Couplings between the horizontal and vertical 
branches of the contour are determined by correlation 
functions and G^. These functions take into ac¬ 
count the coupling “across the quench”. Notice, that the 
Matsubara field ipM only couples to the quantum com¬ 
ponent For the bare Green’s functions this implies 
that only the Keldysh function contains pre-quench mem¬ 
ory effects. Thus, the bare post-quench Keldysh function 
g^ {t, t') of a non-interacting system will depend on both 
time variables, i.e., display aging effects. How to deter¬ 
mine g^ {t, t') will be discussed in detail below. The bare 
retarded function only depends on the time differences 


We consider interaction effects within the large-iV ap¬ 
proach. For a large number of field components N 
the generating functional W\h] can be evaluated in the 
saddle-point approximation controlled in small 1/iV. In 
Appendix B, we explicitly perform this approximation 
on the three-time-contour and summarize the main steps 
of the large-A^ analysis for a system with quench. The 
resulting self-consistent large-A^ equations before the 
quench read 

hi = ri(j)i (41) 

Ti = ro,i -b +Ui [ Gf. {q,ujn) (42) 

^ J q,LOri 

and the equations after the quench read 

noo 

hf = — dt'{G^)~^ 

Jo 

—4>i dt'Sr] {t — t') (43) 

J —oo 

(i) = ro./ + ^ id, t, t). (44) 

We use the notation f = f td and f = T V for the 
momentum integration and the summation over Matsub¬ 
ara frequencies, (jii and are the order parameter and 
renormalized mass prior to the quench. G^ is the renor¬ 
malized Matsubara function, where fo,i in Eq- (31) is re¬ 
placed by ri- In full analogy is the renormalized retarded 
function 

+ Sr]{t- t') (45) 

governed by the renormalized, time-dependent mass 
r{t). The non-equilibrium dynamics leads to a time- 
dependence of the mass such that the response function 
is affected by aging behavior. The set of equations is 
closed by the Keldysh function that can be expressed in 
the form 

G^{t,t')= [ G^it,s)M{s,s')G^{s,t'), (46) 

J s,s' 




where we introduced the memory function M (s, s') that 
is explicitly determined in Appendix C. Once this mem¬ 
ory function is known we have a closed set of equations for 
the non-equilibrium dynamics after a quantum quench. 

For our subsequent analysis it is convenient to rewrite 
the last two equations in terms of an expansion in 


5r {t) = r (t) — rf , (47) 


where rf = r (t ^ oo) is the equilibrium value of the 
renormalized mass for the post quench parameters (see 
Eq. (53) below). This includes a shift due to interactions. 
It implies in particular that r/ = 0 for a quench right 
to the critical point. The renormalized propagators can 
then be obtained from the Dyson equations 

(48) 

-I- dsg^ (t - s) Sr (s) Gf (s, t') 

and 

= g^ (49) 

-I- f dsg^ {t, s) Sr (s) G^ {s,t') 

Jo 

+ f dsg^ {t — s) Sr {s) G^ {s,t').. 

Jo 

Here, the function g^ {t — t') is the solution of Eq. (45) 
with r (t) replaced hy rf, i.e. it reads in Fourier space as 




1 

Lu^ — rf — + Sr]{uj) 


(50) 


The Keldysh function g^ is determined as in Eq. (46) 
via g^ {t, t') = g^(t - s)M{s, s')g^{s' - t'). 

In Secs. V and VI we present solutions of the large- 
N equations for several different quench protocols, both 
starting in the disorder'd and in the ordered phase. In 
the next sections, however, we first summarize the equi¬ 
librium limit without quench and motivate the scaling 
behavior out-of-equilibrium of various observables that 
will be confirmed by the subsequent explicit calculation. 
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it follows with dtSrj (t) = 0 that the equations be¬ 
fore and after the quench give, as expected, the same 
solution. We can drop the subscripts i and / that dis¬ 
tinguish between the pre- and post-quench regimes. The 
renormalized mass fulfills 

r = fo + uf G^{q,uj)+ ^ 41 ^ . (53) 

The equation of state for the order parameter reads 

r4> = h. (54) 


In the ordered phase it holds that 0 7 ^ 0 and using 
Eq. (54) it thus follows for = 0 that r {h = 0) = 0. 
The excitation spectrum is massless due to the Gold- 
stone theorem and the fact that longitudinal excitations 
are of order 1/A^. The transition takes place when r = 0 
and 0 = 0 and determines the critical value fo,c (u). For 
T = 0 the latter is determined by 


ro,c = -u 



1 

0}“^ + q^ + Sr]^ {(jj) ’ 


(55) 


depends on the momentum cut off A and, via rj (0), 
on the frequency cut off of the bath spectrum. However, 
universal behavior emerges as a function of the distance 
Sr = fo — ro,c to the quantum critical point. For suffi¬ 
ciently small energies the dynamics is dominated by the 
bath. Comparing length and time scales in the prop¬ 
agator with Sr]^(uj) oc \uj\°‘ yields the dynamic scaling 
exponent 


z = 2/a. (56) 

Analyzing Eqs. (53) and (54) one obtains below the upper 
critical dimension d < due = 4 — z the well known results 
for the order parameter |0)| oc {—Sr)^ with /3 = 1/2, for 
the correlation length ^ oc Sr~'' with l'~^ = 

d + z — 2 , and for the order-parameter susceptibility at 
the critical point G^ (q = 0 , w = 0 ) oc Sr~'^ with 7 = 2iy. 
For d > due all exponents take their mean field values. 


IV. SCALING BEHAVIOR 


C. The equilibrium limit of the large-N analysis 


Without quench, the system remains in equilibrium 
and the large-fV equations before and after the quench 
become identical. In equilibrium, the Keldysh function 
G^[q,t,t') only depends on the time difference t — t', 
with Fourier transform determined by the fluctuation- 
dissipation theorem: 

Gf (g,w) = 2fcoth(^^)lmG/^((7,a;). (51) 


Using 


tJ2g^ {q, w„) = i y ( 9 , 


(52) 


In this section we discuss scaling arguments for the 
non-equilibrium dynamics after a quench towards the 
critical point. In addition to the dynamics of the or¬ 
der parameter 0 (t)</’ = {^p {k, t)) we are interested in the 
retarded and Keldysh Green’s functions G^ (fc, t, t') and 
G^ (k, t, t'), respectively, where both time arguments are 
after the quench. We suppress the field index I here and 
in the following for simplicity. These scaling relations and 
the corresponding exponents will be determined from a 
self-consistent solution of the coupled large-V equation 
in the subsequent sections. 

In equilibrium, the order-parameter at T = 0 obeys 

0eq (Sr, h) = 6-^/"0eq , (57) 
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with distance to the critical point Sr = tq — fo,c and 
scaling parameter b. Note, Sr can in fact be realized by 
either changing fg or by changing fg.c that depends on u. 
Thus, for the discussion of scaling behavior, it is sufficient 
to include only two scaling fields R = {Sr, h). The choice 
b = Sr~'^ implies for h = 0 that (peq (Sr) oc Sr^, where 
/3 is the order-parameter exponent, v is the correlation 
length exponent with equilibrium correlation length 

^ oc Sr~‘'. (58) 

In analogy, right at the critical point and for finite field 
follows with b = that (/)eq (h) oc . Those are 

the well known scaling relations in equilibrium. 

In a non-equilibrium setting the order parameter will, 
on the one hand, depend on time, which transforms un¬ 
der scaling as t —^ b~^t. Here, 0 is the dynamic scaling 
exponent that relates typical time and length scales. We 
consider a regime where the dynamics is dominated by 
the coupling to the bath. In this case, to leading order in 
\/N we find that like in the equilibrium case discussed 
previously (see Eq. (56)) that 

z = 2/a (59) 

is determined by the spectral function rj (w). On the other 
hand, our quench protocol implies that the order param¬ 
eter depends on the initial distance Sri and the final dis¬ 
tance Srf to the critical point as well as on the initial 
and final fields hi and hf. It just holds that 

{t, R, (1), Rf (1)) = 6-^/^'/. {b-H, R, (6), Rf (b)) , 

(60) 

where Rf {b) = (b^^''Srf,b^^^'^hf) and Ri{b) = 

{y^l''Sri,b'^^^l''h^. The scaling dimension of the order 
parameter and of the final values Srf and hf continue 
to take their equilibrium values (3, and fSS/u respec¬ 
tively. This is a consequence of the fact that the system 
approaches equilibrium after very long time scales 


with b~^t = and tj being a microscopic time scale. 
We can therefore express the order parameter in terms 
of a scaling function as 


(j){t,Sri) = t ^‘'^/{t/t*), (63) 


where t* cx Sr^ and t, t* are in units of tj. For t ^ t* 
the scaling function is expected to approach a constant 
'I'(a; » 1 ) —>• const, and the order parameter decays adia- 
batically according to a power-law characterized by equi¬ 
librium exponents 

(64) 


On the other hand, for short times t ^ t* one expects 
that (j) oc (5rf such that 'k(a: 1) oc In this 

regime follows 


</(t<r)oct("-i)^/(""). (65) 


As shown below the exponent k varies with z and is for 
some z-values larger than unity, implying that the order 
parameter rises for intermediate time scales, while for 
other z values k < 1 such that a slower decay of gov¬ 
erns the regime up to the crossover scale t*. Note that 
t* diverges for shallow quenches where Sri is small. An 
analogous behavior occurs in the regime where we per¬ 
form a quench of the field hi to hf = 0 for a system at 
the critical point: Sri = Srf = 0. The characteristic time 

scale is now given as t* cx h^ . The resulting time 
dependence for t ^ t*and t t* are the same as before. 

Let us now turn to the correlation and response func¬ 
tions. In equilibrium, the retarded and Keldysh Green’s 
functions only depend on the difference t — t' ot the two 
time variables and obey at the critical point the estab¬ 
lished scaling behavior 


{q, t, t') = 




R{K) 


eq (g- 

g2-z-r)^z/2 


( 66 ) 


lim (/(t,i?i,i?/) = (/eq(^/), ( 61 ) 

t—¥OQ 

independent on the initial values Sri and hi. Thus, the 
corresponding scaling dimension are the same as in equi¬ 
librium. There is, however, no reason why the scaling 
dimension of Sri and hi should also be equal to the equi¬ 
librium values. This is reflected in the new exponent k 
that enters Eq. (60). At this point it is not obvious why 
the same exponent k modifies the scaling dimension of 
Sri and hi. This will only become clear in our explicit 
analysis of the theory. To demonstrate this and to deter¬ 
mine K is one of the goals of this paper. 

Let us further analyze implications of this scaling be¬ 
havior. For simplicity, we first consider the case of a 
quench right to the critical point Srf = 0 in vanishing 
field hi = hf = 0. In this case Eq. (60) simplifies to 

4> {t, Sr,) = b-^/^cS {b-% b^/'^Sr^ . (62) 


with scaling functions {x) and {x). The factor 
occurs to render the scaling functions and their ar¬ 
guments dimensionless. As discussed, this behavior is 
only valid for time scales sufficiently long to ensure that 
the dynamics is dominated by the bath and it is justified 
to neglect compared to |cli|“ in the bare propagators. 

Considering now the out-of-equilibrium dynamics after 
the quench to the critical point, correlations and response 
depend in general on both time variables. This gives rise 
to the additional dimensionless ratio t/t', compared to 
scaling in equilibrium. Thus, we expect from dimensional 
considerations the following behavior: 


{k,t,t') 


/ t \ {qHl^G, t/f) 

\t'/ g2-z^z/2 


Here, the scaling functions F^ and F^ are chosen such 
that they depend only weakly on the ratio t'/t if t ^ 
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t'. The possibility of a singular dependence on t/t' is 
captured by the exponents 9 and 9' ^ respectively. Below 
we will determine scaling laws that relate the exponents 
0 , 0 ', and K in such a way that only one new independent 
exponent emerges. Thus, the post quench dynamics is 
governed by a single critical exponent that cannot be 
expressed in terms of equilibrium exponents. 


V. QUENCH FROM THE DISORDERED PHASE 

In what follows we present the solution of the cou¬ 
pled large-A equations (see Eqs. (41)-(46)) for the post- 
quench non-equilibrium dynamics at the quantum criti¬ 
cal point starting the quench from the disordered phase. 
The quench protocol is indicated in by A —>■ C in Fig. 2. 
The effective time dependent mass in Eq. (44) simplifies 
to 

r{t) = roj + ^ [ (g,t,t). (68) 

^ J (27r) 

In case of a quench to the critical point, tqj takes the 
value fo,c in Ed- (55). It is convenient to express fo,c in 
terms of the equilibrium Keldysh function such that 

^ J (27r) 

For times larger than the microscopic time scales the 
Keldysh function obeys scaling. We refer to Appendix A 
for a discussion of the various short-time scales in the 
problem. Based on dimensional arguments, we make the 
following generalized light-cone ansatz 

r(() = ^. (70) 

valid for t larger than the microscopic time scale. The 
prefactor 7 , that determines the strength of the coupling 
to the bath, was chosen to ensure that the light-cone 
amplitude a is dimensionless. This ansatz corresponds 
to a time-dependent correlation length 

^(^) = |r(^)r^ = ^^^^l/^ (71) 

(7|a|)i/2 

reflecting the fact that after time t perturbations propa¬ 
gate a distance ^ (t). For ballistic propagation with z = 1 
this corresponds to the usual light-cone propagation dur¬ 
ing time t. In general, ^(t) is the length scale on which 
the signal of a perturbation can propagate. Among oth¬ 
ers, this implies that the system will eventually reach a 
state with diverging correlation length. This is expected 
as the coupling to the heat bath ensures that the system 
equilibrates at the quantum critical point since the bath 
temperature is zero. 

We now demonstrate that Eq. (70) indeed leads to a 
self-consistent solution of the coupled large-A equations 
and determine the light-cone amplitude a. Let us next 
explore some consequences of such a time dependent mass 
r(t). 


A. Light-cone amplitude a and universal 
short-time exponent 6 


We first establish a relation between the light-cone am¬ 
plitude a of Eq. (70) and the exponents 9 and 9' of the 
retarded and Keldysh Green’s functions in Eq. (67). We 
will show that 9 = 9', which is in contrast to the regime 
near classical phase transitions^^’®^ and to the dynamics 
in an isolated, quantum system®®. This relation can eas¬ 
ily be determined from an analysis of the Dyson equation 
of in the intermediate regime tj t tq = (i/y/g) 
and t ^ t'. The leading corrections that follow from 
Eqs. (48) and (49) in this limit are given by 

6G^ {k, t, t') « dsg^ {k, t - s) {k, s - t'\12) 

SG^ (k, t, t') « ^ dsg^ {k, t - s) ^5^ {k, s, t') (73) 

Here, the lower integratl limit appears since the power- 
law decay of r {t) oc sets in only beyond this time 

scale. 

We first consider the correction to the retarded func¬ 
tion. For frequencies uj >> {q/s/y) , which correspond to 
times t <C , the bare retarded propagator g^ at 

the critical point is local in space and entirely governed 
by the bath 


9’^{q,oj) 


I 

5iq{uj) ■ 


(74) 


Here, 5rj (w) is the frequency dependent contribution of 
the bath function rj{uj) (see Eq. (12)). Performing the 
inverse Laplace transformation we obtain 


9^ iq,t) 


sm{Tr/z) 2 /z-i 
7 r(2/z) 


(75) 


If we insert this into Eq. (73), it follows for t' > that 


SG^ 


asin^(7r/z) f* — ^ {s — ^ 

'yT{2/zY Jt, ^ s2/z 


The integral can be performed analytically. Evaluating 
it in the limit t ^ t' yields 


SG^ 


asin^(7^ 2/z-i 

7r(2/z)2 



(76) 


To leading order it follows for the retarded propagator 
G^ {t, t') = g^ {t — t') ^1 -f 0 log — -I- .. .^ , (77) 

with 


a sin(7r/z) 

r(2/z) • 


This result demonstrates that a time-dependent mass 
r{t) leads to aging behavior of the retarded Green’s func¬ 
tion, i.e. G^{t,t') depends on both time arguments t 
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and t' separately. If we include higher order corrections, 
an analysis along the same lines yields higher powers of 
the logarithm with appropriate coefficients allowing us 
to exponentiate the logarithm if a is small. We can then 
put {t, t') in the scaling form of Eq. (67) and find 
that 9 is indeed the exponent that appears in the scal¬ 
ing form. Thus, the dimensionless amplitude a of the 
time-dependent mass determines the exponent 9 of the 
retarded Green’s function. This analysis was performed 
under the assumption of negligible g—dependence of the 
retarded propagator. Including this q—dependence is 
technically slightly more involved, but leads to the same 
behavior since the log(t/t')-time dependence is still the 
dominant one in the limit 

Let us now discuss the interaction correction to 
the short-time behavior of the Keldysh function using 
Eq. (73). We need to evaluate the short time behavior of 
the non-interacting Keldysh function for which 

it is useful to introduce a double Laplace-transformation 
(see Ref. 81) 

noo 

Jo 

(79) 

As we show in Appendix C the Keldysh function can be 
expressed in terms of the memory function M[q,ui^uj') 
as g^{q,u},uj') = M{q,uj,uj')g^{q,uj)g^{q,uj'). Using the 
deep quench result for the memory funcion in Eq. (C16) 
at r = 0, this takes the form 

Kf _ .sign(w)d 77 (w) -k sign{u}')6r]{uj') 

X ff^(9,w)5^(g,a;'). (80) 

The limit t ^ t' corresponds to w' ^ w. In this limit 
5rj(ui') is always large compared to 6 g{u}) and it follows 

g^{q,UJ,u;' » w) ~ \ ^{q,uj)g^{q,uj'). (81) 

Iw I 

Using Eq. (74) for the retarded function in the short time 
limit, it immediately follows 

g^{q, uj, w' » w) ~ Y^^g^{q, uj). (82) 

|w'| 


B. Long-time behavior of 

We next determine the long time decay of the Green’s 
functions and G^ of a system with time-dependent 
effective mass r{t) as given in Eq. (70). To this end we 
analyze the equation of motion (43), in the long time 
limit. From the equation of motion follows 

{dl +r{t) + q^) ip {q, t) = f dsSg{t - s)p{q, s). 

J —OO 

(84) 

In order to analyze the impact of the time dependence of 
the correlation length, we perform an expansion for small 
r (t). For the specific ansatz in Eq. (70), this amounts to 
an expansion in the dimensionless light-cone amplitude 
a: 

¥’(9,0 = ¥’eq(9,0 + ¥’ 1 ( 9,0 + C’(a0, (85) 

where ¥’eq(9,0 obeys the equation of motion for r (t) = 0, 
i.e. the equation of motion in equilibrium at the critical 
point. We can split the integral on the right hand side of 
Eq. (84) into 

H (t) = f ds5f]{t — s)p{q, s) (86) 

J —00 

and a contribution that only contains the post-quench 
dynamics of p {q, t). As far as the post-quench dynamics 
is concerned, H {t) acts as an inhomogeneity in Eq. (84) . 
At long times we can neglect this inhomogeneity since 
Sr] {t ^ 00 ) vanishes sufficiently rapidly. This determines 
the leading perturbative correction to 

¥’ 1 ( 9,0 = -/' dsG^qit-s)r{s)p^qiq,s), (87) 

Jto 

where Gfq(t-t') = -i9{t-t'){[peqit), peqit')]) (suppress¬ 
ing the field index I again). 

Using the definition of G^^k |' 3 g^ .^g 

can now determine how the Green’s functions relax to¬ 
wards their equilibrium expressions. Let us first consider 
the equal-time Keldysh function G^{q,t,t), which turns 
out to be of particular importance for the self-consistent 
solution of the large-A^ equation. Using Eq. (85) leads to 


The back transformation is now straightforward and, us¬ 
ing that l/|w'| turns into a constant upon Laplace trans¬ 
formation, we find 

5 ^(g,t » t') c. g^{q,t) = (83) 

7r(2/z) 

In the short time limit and for t ^ t', the Dyson equa¬ 
tions for the retarded and Keldysh functions thus only 
differ in the bare values of the two functions. Since 
we just demonstrated that g^{q,t t') and g^{q,t) 
have the same t—dependence in the relevant regime, both 
Dyson equations have the same solutions. This immedi¬ 
ately yields 9' = 9 for the exponent 9' of the Keldysh 
function introduced in Eq. (67). 


G^iq,t,t) = G^q{q) + 6G^ {q,t,t), (88) 

where 

SG^ {q,t,t) =-2 f dsG^^it-s)r{s)G^^it-s). (89) 

To evaluate this expression we express the retarded func¬ 
tion in Laplace and the equilibrium Keldysh function in 
Fourier space and obtain 

5G^ {q,t,t) = -2i J ^^^^Im[G^(a;)]G^(a;') 

nt 

X / dsr(s)e-*(“-“')(‘-") . (90) 
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We want to evaluate this integral in the limit of large 
t. The typical frequency scale of the Green’s function is 
^ For times t ^ the integrand 

is in general highly oscillatory, except for contributions 
that stem from the upper limit of the integration over 
time. This yields up to leading order in t 

r (s) « . (91) 

Jt^ w - w' 

Performing the remaining frequency integration leads to 


Gr iqG,t) = Geq(g) + 


2 ir (t) 




with numerical coefficient 


ck = 


4sin(7rz/2) 


z{2 — z) sin^/^(7r/z) 


(92) 


(93) 


This result demonstrates that the equal-time correlation 
function decays for large times according to a power 
law. Including higher order corrections in the Dyson 
equation, leads to terms of order r{t)^. Thus, the term 
= 2ir{t)f(^CKq‘^~^J^^'^) appearing in Eq. (92) 
is the slowest decaying correction at large times. This is 
an interaction effect, because the bare correlation func¬ 
tion g^{q,t,t) decays for finite momenta always expo¬ 
nentially oc exp to its value in equilibrium. 

Critical fluctuations lead to a significant slowing down of 
the equilibration of the system. 


C. Self-consistent determination of the light-cone 
amplitude 


decay much more slowly. Using Eq. (92) such corrections 
G^{q,t,t) oc r{t)lq'^~^ are proportional to Upon 

integration this generates terms that behave as 1/e for 
small e: 





(96) 


where qo is some appropriate lower cut-off that we elab¬ 
orate on below. Such a term, multiplied with r (t) oc a 
will be of same order 0(e°) as the bare . Thus, in 
the expansion of G^ in a we have to keep those slowly 
decaying terms proportional to . Erom Gf in 

Eq. (92) follows that large times t correspond to large 
momenta q ^ qo = ^ which naturally intro¬ 

duces the lower cut-off qo in Eq. (96). As we have pre¬ 
viously shown in Sec. V A, in the opposite limit of small 
momenta q the Keldysh-function becomes 

momentum independent (see Eq. (83)) and cannot gen¬ 
erate contributions that behave as 1/e. 

With these considerations, we are able to solve the 
large-A^ equations self-consistently. We expand G^ = 
g^ + G^ in small r{t) to find 


r{t) 




uKdr{t)/ 

67-/2 V V fi/- J ) ■ 


(97) 


In the first integral we already took the limit A —>■ 00 be¬ 
cause g^ approaches G^ exponentially quickly and the 
integral is thus convergent at the upper limit. To pro¬ 
ceed, we use the scaling form of Eq. (66) and introduce a 
dimensionless integration variable x = q^tj^^G to obtain 


In this section, we demonstrate that the ansatz for the 
time-dependent mass r(t) in Eq. (70) indeed leads to a 
self-consistent solution of the coupled large-equations 
and we show how to determine a. Let us begin with 
some general remarks about Eq. (69) that we want to 
solve self-consistently 

^ = ^^jo dqq^-^-^ {iG^(q,t,t)-iG^^{q)) . 

(94) 

Here, Kd = r(d/2)/[27r‘^/^(27r)'^] takes into account the 
integration over angles and 

e = 4 - d - z . (95) 

We want to determine the dimensionless light cone am¬ 
plitude a self-consistently for small e. Since u will turn 
out be also of order e (see Eq. (100) below) it is tempt¬ 
ing to speculate that one only has to expand G^ in a. 
To leading order it would then suffice to insert the bare 
Keldysh function g^{q,t,t) in Eq. (94). Since g^ de¬ 
cays exponentially to G^ the integral is convergent in 
the limit A —>■ cx). However, as was shown in the previous 
paragraph, first order corrections to the Keldysh function 


«7 _ uKd jGotG^ ua'yKd / _ 

t^/- 2 z 7-/^ 7''/-f^/- ’ 

(98) 

where we have introduced the dimensionless integral 


pOO 

Go = / (i/^(a:, 1). (99) 

Jo 


The equation (98) must hold for all times t ^ tj, and 
therefore terms with and without U/- must cancel sepa¬ 
rately. This allows to extract the two conditions 


* _ ckT 
u = u = 


72 A6 


Kd 


CkGq 

a = — —e. 
2z 


( 100 ) 

( 101 ) 


The first equation determines the deep-quench fixed 
point value of u. It does not imply that u must be tuned 
precisely to this value. Instead, it should be understood 
as a consequence of using the deep-quench scaling form 
of G^. Non-universal contributions in G^ will renor¬ 
malize It, which is fully consistent with our renormaliza¬ 
tion group based reasoning in Ref. 82. The second equa¬ 
tion (101) determines the value of a and, using Eq. (78), 
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d je 



Z 

Figure 4. (Color online) Prethermalization exponent 0 as a 
function of dynamic critical exponents a. Here, e = 4 — d — 2 
and 2 = 2/a is determined by the form of the bath spectral 
function at low energies Imr 7 (a;) oc |a;|“. For 6 < 0 one ob¬ 
serves an underdamped approach of the free Keldysh Green’s 
function towards equilibrium. The transition from 

overdamped to underdamped behavior occurs close the loca¬ 
tion of the sign change. 


the value of the new critical exponent 


sin(7r/2)a 

r(2/z) 


CK sin(7r/2) 

2zr{2/z) 


( 102 ) 


Note, that exactly the same value was obtained in our 
previous renormalization-group analysis in the limit N —>■ 
oo®^, where we have shown that the integral Cq can be 
performed analytically in the case of an ohmic bath. For 
general dynamic critical exponents 1 < 2 < 4 we evaluate 
it numerically and the result is presented in Fig. 4. 

The relation between the light-cone amplitude a and 
the retarded short-time critical exponent 9 is identical 
for a classical and a quantum system, however, the 2 - 
dependence of 9 is completely different due to the fact 
that quantum fluctuations are only present in the quan¬ 
tum system. Quantum fluctuations influence 9 in two 
ways: (i) via the numerical constant ck in Eq. (93) and 
(ii) via the result of the dimensionless integral C'o over 
the Keldysh functions in Eq. (99). 

The constant Ck determines the fixed point u* as well 
as the long-time decay of . For a quantum system, 
Ck is positive for 2 < 4 and linearly approaches zero for 
2 —4. The constant Cq is given by the difference of 
the post quench Keldysh scaling function f^{x) and the 
equilibrium Keldysh scaling function F^. For a classi¬ 
cal system, the system is always overdamped and is 
smaller than for all values of 2 . Hence, the expo¬ 
nent for a classical system is always positive 9ci > 0. For 
an ohmic bath, for example, one finds 9ci = Cci/4with 
Cc/ = 4 — For a quantum system, however, can 
exhibit oscillations, in particular for super-Ohmic bath 


spectral functions where 2 < 2 the Keldysh function be¬ 
comes underdamped. This leads to a sign change of 9 at 

2 « 1.8. 

While the classical result of the exponent for an ohmic 
bath is 9ci = which is the same value as in the 

quantum case for 2 = 2: 9{z = 2) = e/4, this is only a 
coincidence, since the constants ck and C'o take different 
values for T = 0 and T —>■ 00 . For an ohmic bath, one can 
therefore still distinguish the quantum post-quench dy¬ 
namics from the classical one by analyzing, for example, 
the long-time behavior of the Keldysh-function. More¬ 
over, since e depends on 2 , we find the same short time 
exponent for classical and quantum systems in different 
dimensions. For example, we predict 9^1 = 1/4 for a clas¬ 
sical system in contact with an Ohmic bath in d = 3 
dimensions, while 9 — 1/4 for a quantum system in con¬ 
tact with an Ohmic bath in d = 1 dimension. 


D. Distribution function n{t,uj) 

The analysis of the long-time behavior of the equal¬ 
time Keldysh function G^{q,t,t) in Sec. VB can be 
straightforwardly extended to different times t, t' as well 
as to retarded function G^{q,t,t'). We show in detail 
in Appendix D that in the limit where both time ar¬ 
guments are large compared to the typical mode time 
tq = but the relative time is small t — t' tq, 

this leads to 

Gf (d, t, t') =G^q{q, t-t')- 2r C^{q, t - t') 

(103) 

G^{q,t,t')=G^^{q,t-t') (104) 

- 4i9{t - t')r G^(g, t - t') 

where we have defined the functions 

=y ^Re[Gfq(g,a;)]G^((7,w)e"“* (105) 

G^(d,t) = J ^Re[Gfq(g,a;)]lm[Gfq(q,a;)]e-“*. 

(106) 

By performing a Wigner-transformation of Eqs. (103) 
and (104) and by using the fluctuation-dissipation the¬ 
orem for the equilibrium Keldysh function G^ = 
2i coth[w/(2T)]lmG^, one can express G// after a few 
steps as 

G^{q,ta,uj) = 2icoth [l -b 2r(ta)ReG^q(q, w)] 

X lmG^(g,ta, w), (107) 

where ta = {t +1')/2. This result shows, that we are in 
the limit of adiabatic relaxation, since G// can be written 
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G^{q,ta,uj) =coth 

2 ilm? 7 (a;) 

+ r{ta) + Re?7(w)]2 + [Iin77(w)]2 ‘ 

(108) 


Further, the similarity between and G^ suggest to 
connect them via the fluctuation-dissipation theorem and 
to introduce a distribution function n(ta,uj) via 

Gf (g, ta,uj) = 2i [2n{ta,uj) + 1] ImG^(g, ta,u}). (109) 


With n(ta,ui) = nsiuj) + 6 n(ta,uj), where nBi'jj) is the 
Bose-distribution function, this yields for the correction 

Sn{ta,uj) = coth r(ta)ReGfq(g,a;). ( 110 ) 


This expression is only valid for large frequencies, which 
corresponds to small relative times \t — t'\ <C tq. For 
ta —t oo the system relaxes to the thermal equilibrium 
state. Aging effects, however, lead to a significant slowing 
down of thermalization. Explicitly, one finds 


6 n = coth 



9Ti2/z) - 

\2T) 

(|cc|G) 2 A^ 


TT 

cos —h 
z 


( 111 ) 


where both q^\uj\ = {\uj\tq) 1 as well as 

< 1 , since ta'> tq'> \t - t'\. 

The correction to the Bose distribution is proportional 
to the short-time exponent 6 and exhibits a (dominant) 
algebraic frequency dependence 6 n oc for z 2 

and 5n oc in the ohmic case. This slow decay at 

large frequencies clearly shows that the system is not 
thermal and cannot be characterized by a temperature 
as the decay would then be exponential in w. The sign of 
5n is positive for a (sub)-Ohmic bath z>2, showing that 
even long after the quench there are an increased number 
of excitations at large frequencies present in the system as 
compared to equilibrium. For super-Ohmic bath spectral 
densities, however, the sign of 5n can change depending 
on the product of 9, which becomes negative for z < 1.8, 
and the last term in Eq. (111). This term stems from 
Re G^ and is positive (negative) for z > 2 (z < 2), except 
slightly below z = 2 where the sign depends on the in¬ 
terplay of the cos(7r/z) < 0 and (|a;|tg)“^A sin( 7 r/z) > 0 
with 1. More importantly, a negative sign 

of 5n implies that the density matrix of the system is not 
diagonal in the energy basis. The off-diagonal terms de¬ 
scribe the presence of quantum coherence and thus make 
a straightforward interpretation of 5n in terms of a dis¬ 
tribution function impossible. 

In Fig. 5 we illustrate this result for 2n{ta,uj) -|- 1 in 
the ohmic case, both for 6* > 0 as well as for 6* < 0 
(although 6 = e/4 > 0 for z = 2). We include 0 < 0 in 
the plot to illustrate the behavior in cases where 9 < 0 
which is qualitatively the same. Above the upper critical 
dimension, one finds 9 = 0 and thus Sn = 0. 


coth(^)(l -b 2r(ta)ReG^J 



Figure 5. Visualization of the correction to the distribu¬ 
tion function 5n{ta,uj) as a function of dimensionless vari¬ 
able uJ'y/q^ (see Eq. (110)). Plot shows 2n{ta,oj) -I- 1 = 
coth(/l/)(l -I- 2r{ta)ReG^q) (see Eq. (107)) as a function of 
the dimensionless variable x = for z = 2, jT/q^ = 1.7 

and | 6 l| 7 /(g^taG) = 1 . 


E. Upper critical dimension 


From the bare Green’s function in Eq. (40) it fol¬ 
lows immediately that the mean-field value of the ex¬ 
ponent 9 is zero. Therefore, 9 vanishes above the upper 
critical dimension due = 4 —z. In this section, we investi¬ 
gate the post-quench behavior if the system is exactly at 
the upper critical dimension d = due- We can straight¬ 
forwardly extend our previous result to this case. The 
large-A^ equation for the mass reads 

Tucit) = j dqq^-^ [fG^^(g,t,t) - fG^(g)] , 

( 112 ) 

where we used that d = due- In a next step we expand 
G^^ in Tucit). The important slowly decaying correction 
comes from terms going with G^{q,t,t) oc q~‘^~^^, gener¬ 
ating logarithmic divergences after the g-integration. In 
Sec. VB we calculated this correction, without making 
any assumption about the dimensionality of the system 
or the concrete decay of the time-dependent mass r(t). 
Therefore, one can directly use the result of Eq. (92). 
Again, the condition t > tq = translates into 

a lower cut-off in the g-integral. For q < q^ = 
the Keldysh function G^^ is again local in space, inde¬ 
pendent of the dimensionality, and cannot generate loga¬ 
rithmic terms. Focusing on the contribution of momenta 
g > go it holds that 


/q,l/2p/z 


dgg'^ ^fGf(g,t,t) 


2 ruc{t) . / 

7I/2 j 


2 r„e(t) 


In 




(113) 
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Compared to the case d < due where e > 0, the slowly de¬ 
caying correction is now logarithmically divergent both 
in the cut-off A and at the lower boundary q^. This 
is not surprising, because the same behavior occurs in 
equilibrium. The ^-integration of higher order terms in 
the expansion in r{t) does not generate logarithmically 
divergent contributions and can therefore be neglected. 
Inserting the expansion = G^ -I- Gf into Eq. (112) 
and keeping only the logarithmically divergent correc¬ 
tions, we obtain 

^ u-fKdCp uKd / f 1\\ 

2z^GH'^Gruc{t) \ \AH^ ) ^) J ' 

(114) 

This equation holds for all times t > ii contributions 
with and without t-dependence cancel separately. This 
give rise to two conditions, one for the time-dependent 
mass ruc{t) and one for the deep-quench fixed point value 
of u: 


VI. QUENCH STARTING IN THE ORDERED 
PHASE 

In this section we will consider a quench with a finite 
initial order parameter < Q) = (pi. A finite order 
parameter can be achieved either by (i) a finite initial 
field hi, which is switched off at t = 0 (path B' ^ C 
in Fig. 2) or by (ii) a quench in the mass-parameter 
with fo,i < f'o.c (path B G in Fig. 2). Both quench 
protocols yield the same exponent 9 for the retarded 
and the Keldysh Green’s functions. Furthermore, the 
exponent 9 determines the short time dynamics of (p(t) 
and the magnitude of the typical cross-over time t* 
from prethermalized dynamics (p oc to relaxation 
p oc 

Starting in the symmetry-broken phase, the post- 
quench I/iV-equations for a quench to the quantum crit¬ 
ical point rf = 0 read 


ruc{t) 


u = u 


* 

UC 


ICrCq 

zckY^'^ 

Kd\n[yG‘^/(AH^)] 


(115) 

(116) 


Compared to the result below due in Eqs. (100), (101), e 
is replaced hy z/ \n{t/t.y) in the effective mass. Since t ^ 
tj, Tue{t) is indeed small, justifying the earlier expansion 
of G^^ in Tueit). This additional logarithm in the self¬ 
energy correction ^G^^ (t, t') at intermediate times t.y 
t' t tq = implies that 


-In 


i 9ue K In 


In ( — 


In ( — 


'‘1 / i 


(117) 


^ ~ (i)) 

( 120 ) 

pi I dt'dr] {t — t') = r{t)p{t) — ( ds5rj{t — s)p{s). 

( 121 ) 

In the last line we dropped the dfp term, since we concen¬ 
trate on times t ^ where the dynamics of the system 
is dominated by the bath. We also assume a spatially 
homogeneous order parameter p{t). 

In analogy to the quench starting in the disordered 
phase, we expand in G^ -|- Gf for small r{t) at long 
times, where Gf is given in Eq. (92). This yields 


u 





z/2 


with 


This equation is only fulfilled for all times t ^ t-y if 


9 


UC 


jckGo sin(7r/z) 


r(2/z) 


(118) 


Since the logarithmic corrections are small, we can ex¬ 
ponentiate 9ue- This yields, instead of an algebraic di¬ 
vergence as in Eq. (67), a logarithmic divergence in the 
scaling function 




1 

g2 — 




(119) 

The same behavior can be extended for the Keldysh- 
function and its scaling form at the upper critical di¬ 
mension. 


Uf =u* (123) 

with u* and a given by Eq. (100) and Eq. (101). The 
first equation determines the deep-quench fixed point of 
u. The value of u* is not affected by the quench direction 
and is idental for A, B, B' —> G. The second equation 
determines the effective, time-dependent mass. The first 
term in r{t) is exactly the same as for the quench start¬ 
ing in the symmetric phase. However, the mass is now 
modified by an additional term proportional to p"^. This 
time-dependent mass can be inserted into the equation of 
motion (121) oi p leading to a non-linear, inhomogeneous 
differential equation. 























16 


In the following two subsection we present a solution 
for in the prethermalized regime and in the long time 
limit. The special case where the system is at the upper 
critical dimension d = due will be treated in section VIC. 


A. Order parameter coarsening in the 
prethermalized regime and cross-over time scale t* 

At intermediate times, one can assume that the order 
parameter </> is small due to the collapse of the correla¬ 
tions after the quench. This collapse follows a posteriori 
from our solutions which also predict a collapse of the 
correlation length to a microscopic lengthscales right af¬ 
ter the quench. For a fast quench, the collapse occurs 
on microscopic timescales. In this regime, according to 
Eq. (124) the dominant contribution to r{t) is given by 



(125) 

where t* denotes the boundary of the prethermalization 
regime and will be defined below. This yields in the equa¬ 
tion of motion ( 121 ) 

dsSr]{t - s)(j){s) = 

(f>i [ dsSr]{t — s). (126) 

J —OO 

One can solve this equation with the retarded Greens 
function introduced above or via Laplace transfor¬ 
mation. Let us solve it here via Laplace-transformation. 
We expand ^(t) = (po + 4>i in small a and obtain after a 
few steps 

Mt) =Mh) (127) 

/ t 

dsg^{t,s)^^(l)o{s-t^). (128) 

The dominant contribution to the integral in the limit 

t ^ is logarithmically divergent, yielding 

fvpoitj) (^1 + 9 log , (129) 

with the same 9 as in Eq. (78). For small 9 the logarithm 
can again be exponentiated such that 

= Mil) ■ ( 130 ) 

For positive 9 the order increases after a quench, while 
for negative 9 it decays. Note that in both cases the 
order first collapses right after the quench on microscopic 
timescales t~f to a non-universal value < pi. The 

two cases are schematically shown in Fig. 6 

In both cases the order parameter varies slowly in time 
compared to the algebraic decay of Therefore, 


Figure 6. (Color online) Time evolution of the order parame¬ 
ter pit) following a quench from an initial state pi ^ 0. The 
two panels show the cases of (a) 9 > 0, and (b) 9 < 0. First, 
the order parameter collapses due to the quench. After a 
microscopic timescale (see also Appendix A), the order pa¬ 
rameter recovers for a while t.y < t < t* according to pit) oc t® 
for 6^ > 0, before it eventually relaxes to zero for t > t* via a 
power law described by equilibrium exponents. For 9 < 0 the 
order parameter always decays in a power law fashion but ex¬ 
hibits two different exponents for t < t* and t > t*. The inset 
in (b) shows the time dependence of the correlation length 
which diverges after an initial collapse in a light-cone fashion 

e(t) oc 


independent of the sign of 9 there exists a time t* where 
the underlying assumption of Eq. (125) that the second 
term in Eq. (124) may be neglected compared to the first 
is no longer valid. This yields the crossover time t* be¬ 
tween prethermalization and quasi-adiabatic relaxation 


r = 


2\a\Kdt. 


z/(2+2z9-e) 




00 (^7) 




(131) 

By performing a quench with a small amplitude such that 
the system was initially located close to the quantum crit¬ 
ical point and ^o( 7 ) small, the duration of the prether¬ 
malization regime can be tuned to large timescales. Note 
that for 2 : > 2 a large damping 7 increases t*, while 
for z < 2 a large damping shortens the prethermalized 
regime. 


B. Quasi-adiabatic relaxation to equilibrium at 
long times t ^ t* 

At long times the system thermalizes due to the ex¬ 
ternal bath. For a quench right to the critical point, we 
expect the order parameter to eventually relax to zero. 
Thus, for t t*, where t* is the crossover time defined in 
Eq. (131), the order parameter pit) varies slowly in time 
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such that one can neglect the time derivatives of as 
well as the convolution term 


dt'Sriit — 



dt'5rj{t — t') 


0. (132) 


The inhomogeneous term in Eq. (84) quickly decays to 
zero at long times, such that the equation of motion (121) 
reduces to 


r{t-»t*)(j){t)=Q. (133) 


A non-trivial solution of this equation requires that r{t ^ 
t*) = 0. Inserting this solution into Eq. (124) yields for 
the dynamics of the order parameter 


(j){t > t*) 


2 aKd 


gq2/z-l+e/2g^ 


^-(2-e)/(2z) 


/ 2702 !' 


^-/3/(i/z) ^ 


(134) 


where we used the large-fV values of the exponents iz 
and p. The time-dependence of (j) in the long-time limit 
is the usual adiabatic decay and is fully determined by 
equilibrium exponents. This result confirms the scaling 
analysis in Sec. IV. Like in the long time behavior of the 
Keldysh function, the light-cone amplitude a enters as a 
universal prefactor. 


C. Dynamics of at d = due 


We can extend the previous discussion to the case of 
a system at the upper critical dimension d = due- For 
e = 0 the self-consistent solution of the large-TV equations 
is given by (cf. Eqs. (115) and (116)) 


u = u 


* 

liC’ 


ruc{t) 


z f 'JCkCo 


2 Kd ) ■ 


(135) 

(136) 


In the adiabatic regime, the dynamics of </> is determined 
by the condition in Eq. (133), i.e. r{t t*) = 0 which 
yields 


V CKZ-f^/^ J 

where v and /3 now take their mean-field values /3mf = 
I'MF = 1/2. As expected, at d = due the long-time quasi- 
adiabatic dynamics is given by the mean-field values of 
the equilibrium exponents. 

To obtain the crossover timescale t *, we again compare 
the size of the two terms in Eq. (136). In the prether- 
malized regime t < t*, the pP contribution in Tuc is much 
smaller than ( 7 C;gC'o)/[t^^^ and can thus be ne¬ 

glected. In this limit, the effective mass ruc{t) takes the 


form of Eq. (115) and it follows in analogy to the calcula¬ 
tion in Sec. VIA for the dynamics of the order parameter 

4>uc{t) « (/o + Quo In ^In ^ 

« (/o (^In pj . (138) 

Like in the Green’s function, the additional logarithm in 
the effective mass leads to a even slower increase (de¬ 
crease) of order parameter after the quench. The typical 
cross-over time t* is now given by 

(139) 

Here, we have neglected the logarithmic contribution 
ln(t*/t.y)^®"‘=. The typical crossover time for e = 0 is 
now given by the mean-field values of r'MP = Pmf = 1/2 
and 0 mf = 0 in the exponent of (j)o{t^). 


VII. CONCLUSIONS 

In this article we have provided a detailed analysis of 
the post-quench dynamics of a (^"^-model coupled to an 
environment, that is suddenly moved from an equilibrium 
state near a quantum critical point towards the quantum 
critical point. We have employed non-equilibrium quan¬ 
tum held theory in the limit of large-TV, where N is the 
number of held components of (p. We considered differ¬ 
ent values of the dynamic exponent z, which at large-V 
is determined by the bath spectral function that we as¬ 
sume to be of ohmic, sub-ohmic or super-ohmic form. We 
investigate quenches to the quantum critical point start¬ 
ing from both the symmetric and the symmetry-broken 
phases. 

After the quench, the system is characterized by 
an effective, time-dependent mass r{t) = 
with universal amplitude a. The retarded and 
Keldysh Green’s functions obey scaling = 

Away from 

equilibrium, and G^ depend on both time arguments 
t, t' separately, implying interaction induced aging be¬ 
havior. The singular part of the time dependence in the 
short time limit is characterized by a new, universal crit¬ 
ical exponent 9, which is independent of the equilibrium 
exponents. We have explicitly shown that one obtains 
the same value for 0, independent of the quench direc¬ 
tion and the initial parameters. In the long-time limit 
the algebraic decay of r{t) leads to a signiheant slowing 
down of thermalization. Here the value of 9 enters a uni¬ 
versal prefactor, while the aging behavior of G^ and G^ 
can be expressed by equilibrium exponents. 

For a quench with a finite initial order parameter t/i, 
we have identified three different time regimes in the 
dynamics of the order parameter: (i) a non-universal 
regime at short times t < (ii) a prethermalized uni¬ 
versal regime at intermediate times < t < t* and (iii) 
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a regime at long-times t > t* that is characterized by 
quasi-adiabatic relaxation to equilibrium. In the prether- 
malized regime the order parameter fulfills ^(t) oc , 
where 0 can be positive or negative, depending on the 
dynamic critical exponent z, which is determined by the 
bath spectral function. In the long-time limit the dy¬ 
namics of (p is given by equilibrium exponents. The typ¬ 
ical timescale t* separating the prethermalization from 
the quasi-adiabatic regime depends explicitly of the ini¬ 
tial value of the order-parameter pi. It diverges in the 
limit of a weak quench pi ^ 0. This permits tuning the 
prethermalization regime to extend to long time-scales, 
however, at the cost of a small value of the order param¬ 
eter. 

We also analyzed the dynamics of the order parameter 
p{t) and of the retarded and correlation function 
at the upper critical dimension d = due = 4 — 2 . Com¬ 
pared to the dynamics below due the power-law grow (or 
decay) of p{t) oc in the prethermalization regime is 
now slowed down by the occurence of an additional log¬ 
arithm pue{t) OC . In the long-time limit, the 

dynamics of p is given by the mean-field values of equilib¬ 
rium exponents. The special [ln(t/t.y)]^"'= behavior arises 
because at the mean-field level the exponent 9 vanishes. 
Above the upper critical dimension d > due, one cannot 
use a scaling ansatz for . We expect that the or¬ 
der parameter remains constant during the prethermal¬ 
ization regime, which is followed by an adiabatic decay 
with mean-field exponents at long times. 

The unique universal dynamical behavior we report is 
ultimately a consquence of the collapse of the correlation 
function ^ immediately after the quench and its recovery 
by critical coarsening according to ^(f) oc Due to 

the fast quench the system quickly falls out of equilibrium 
and is left in a highly exited state after the quench. Re¬ 
laxation occurs due to interaction and via coupling to the 
external bath, which is we assume to remain at zero tem¬ 
perature. We observe substantial influences of quantum 
fluctuations and quantum coarsening on the post-quench 
dynamics and relaxation. Quantum fluctuations are es¬ 
pecially important for memory effects that the system 
exhibits which lead to a completely different z depen¬ 
dence of the exponent 9 and of the aging prefactor of 
G^ compared to the classical result. Quantum fluctu¬ 
ations can even result in negative exponents 9 < 0 for 
super-ohmic baths. 

In the long time limit, we show that it is possible to 
connect the response function to the correlation func¬ 
tion, like in the equilibrium fluctuation-dissipation theo¬ 
rem. This allows to introduce a non-thermal and time- 
dependent distribution function n{ta,uj) which exhibits 
aging effects, i. e. it depends on the total time ta passed 
since the quench. At a given time, we find n{ta,uj) = 
nsico) + Sn{ta,u}) with 5n oc 0 (|w|to)“^/^[cos( 7 r/ 2 ) -f 
sin( 7 r/ 2 )]. The correction 6 n exhibits an al¬ 
gebraic frequency dependence at large frequencies which 
clearly distinguishes it from a thermal distribution, which 
shows an exponential behavior for T > uj. The amplitude 


of the correction Sn is proportional to the short-time ex¬ 
ponent 9. The sign of 6 n can either be positive or neg¬ 
ative: while 6 n > 0 can be interpreted as an increased 
number of excitations being present compared to equilib¬ 
rium, the fact that 6 n can become negative implies that 
the density matrix of the system is not diagonal in the 
energy basis. The presence of quantum coherence thus 
defies an interpretation of n{ta,uj) in terms of a distribu¬ 
tion function. 

The Young Investigator Group of P.P.O. received fi¬ 
nancial support from the “Concept for the Future” of 
the KIT within the framework of the German Excellence 
Initiative. 


Appendix A: Short-time scales 

Our problem is governed by several short-time scales 
determined by the momentum cut of A of the (p^— theory, 
the strength of the damping coefficient 7 and the upper 
cutoff of the bath spectrum oje- Here we briefly discuss 
the hierarchy of these scales. 

The large-fV analysis clearly shows that we will only 
have a solution r (t) oc if to with 

to = 7^/VAA (A1) 

Otherwise the long time expansion of the upper cut off 
of the momentum integration is not allowed, to obeys 
A^ = 7 t(C^^^, i.e. it is the scale where the damping 
is comparable to the largest possible value of with 
momentum k. Clearly, this comparison is only sensible if 
in fact to > tc with 

tc = 1/wc (A2) 

determined by the bath cut-off. For time scales smaller 
than tc, the damping is very small and it cannot be es¬ 
timated any longer via . The condition to > tc 

translates into 

7^2/" > A2, (A3) 

i.e. even the largest possible value is still smaller than 
the largest possible damping term in the propagator. 

Another assumption of our analysis is that we ignore 
the uf^ term relative to the damping term. This is only 
correct for sufficiently low energies or, equivalently, long 
time scales t ^ t.y with 

= ^-^/[2(^-i)] (A4) 

This result follows from For 2 > I dis¬ 

cussed here, tj can be made small if the damping coeffi¬ 
cient takes large values. However, comparing the ballistic 
and damping term makes only sense if > te, which cor¬ 
responds to 

< Lzl- (A5) 
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The conditions Eqs. (A3) and (A5) imply that it must 
hold 


Wc > A. (A6) 

One can also introduce the time scale 

tA = 1/A (A7) 

which is always larger than tc- Thus, the shortest scale is 
always tc- 

Our prethermalized regime can only start after 
max (to, ^ 7 , ^a)- One finds that no matter what the ratio 
of to and t-y, t\ is always between these two time scales 
and can therefore never be the largest. As for the ratio 
between to and it seems plausible to request ty ^ Ia 
since ty is the only scale that doesn’t depend on a cut-off. 
In this case we have the following hierarchy of scales 

We » A77"/2 » a » (A8 ) 

and universal prethermalized behavior will set in after ty. 
Thus, an important condition for our analysis is that the 
damping due to the bath is sufficiently large such that 
there is a wide regime t ^ ty, where damping dominates 
the dynamics of the system. 


Appendix B: Large-A expansion after a quantum 
quench 

In this appendix we summarize the derivation of the 
large-A equations for the out-of-equilibrium dynamics af¬ 
ter a sudden quench. After having integrated out the 
bath degrees of freedom the action of the collective or¬ 
der parameter field in Eq. (26) is given in terms of the 
Matsubara, quantum and classical field variables as 

S [¥>] = Sm [<Pm] + Sk [<Pc, <Pq] + Sc [<^M, V>q] ■ (BI) 

The first term is up to a trivial factor i the usual Mat¬ 
subara action in equilibrium 

Sm [‘Pm] = l [ Pm (r) (5“)-' (r - r') 9:. (mt') 

+ j (Pm (t) ■ Pm iT-)f , (B2) 

where (r — r') follows from Eq. (31) after Fourier 
transformation. We use the notation J/ = dr with 
usual imaginary time r defined via t = it. For the time 
integrals along the real part we use dt. The two 

segments on the real axis are given as: 





with matrix propagator: 


00 


/- 

9 f 



The inverse is given as 


g 


-1 

0 


0 (ff/)-' \ 

(gf)-i Mo ) ’ 


(B4) 


(B5) 


where Mq = (0o = “(ff/) Since we 

integrated out the bath degrees of freedom the retarded 
Green’s function is determined by Eq. (40) while 


Mq [t, t') = —V (t — t') . (B6) 


The fluctuation-dissipation theorem of the bath can be 
expressed by the Fourier-transforms 

:/(w) = coth(^) (7?(a;)-?7*(w)). (B7) 

The coupling between the two parts of the contour is 


Sc [pm, <fiq\ (t) fj (zT, t) pg (t), (B 8 ) 

J t,T 

where fj {t) = ( 77 ^ (t) — {—t)) and = 4 (jz ± ry =p 

V*)- 

Since ryo = ?y (w = 0) is finite, we have rj (t — t') = 
r]o5{t — t') + 6r] {t — t'), where Sf] [t — t') is purely non¬ 
local. This yields with foj = tqj — rjo that 

(t, t') = - [dt + roj - V^) 6 {t- t') 

+ 5r]{t- t'). (B9) 

Since v (0) = 0, v {t — t') is purely non-local in time. 
Using ry (w = 0) = rjM {icin = 0) we can absorb the same 
zero frequency part on the imaginary axis: fo^i = foy— 70 , 
such that 

(u t') = {pi - fo,i + vp 5 {t- t') 

+ SgM {t - t') . (BIO) 

Next we perform the large-A expansion and follow 
closely the procedure in chapter 30 of Ref. 103. Impor¬ 
tant degrees of freedom are obviously p± = ppl- and 
Pm = ^p'jvi which yields in the case of the round trip 
segment 

1 ^/2 

Pq= ^{P+- P-) = ^Pc- Pq- (Bll) 

The interaction term of the action along the real axis 
part of the contour is then given by 

Skova = - ^ y ^ (Vc (i) + pI P)) ifc {t) ■ pq {t) 

= Pc{t)Pq{t). (B12) 

^ J X,t 
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In order to enforce the above constraints, we introduce 

1= f Dpll6(^{ipl + ipl)-p, 


V2 


IlM IIIM 


x,t 


N‘ 




X e 


-h L,r ^M(t){vl,-NpM) 


(B13) 


where / Dp is the short hand for J Dp^ f Dpq f DpM 
and similar for other degrees of freedom. The integration 
contour for the rc.^-integrations is along the real axis, 
while the rM-integral is performed along the imaginary 
axis. It follows for the generating functional 


(hi^m, 0,..., 0) and the — I components TVc^q^M orthog¬ 
onal to it. 


Vc,q,M — '^c,( 3 ,M^ • 


(BIS) 


We first integrate over the transverse fields ttm along the 
imaginary contour and obtain 

S = Sk + [ (/-M (t) {G^y^ (r, r) (t>M (r) 

+ iN <j)M (r) fj (ir, t) (j)q (t) 

J X,Tt 


+ ^(7V-I)trlog(Gf)-i-^ 


xyt' 


(t') 


-f i- 


.N 


{yPM (t) + fopPM (t) - tm (t) Pm (r)) . 


Z = 




with unchanged coupling term Sc and 


where 

Pit,t')=[ fj{t,iT)Gy {T,T')fj{iT',t') . (BI9) 

J tt' 


Sm = ^ 


V>m{t-){G^) ^ (t,t')<Pm (t') 

't' 

+ *y/ (^PM('r) + (^0x-?’M(T))pM(T)) 


The term f iZqpTVg can be absorbed into a redefinition of 
the bare inverse Keldysh function 

Mo (t, t') M {t, t') = -v{t- t') - p (t, t'). (B20) 


as well as 
1 


Sk = 


! x,ty' 


V’c (t) 




it)) 


N 

~2 


(V2 (ro,/ - rc it)) Pq it) + UfPq it) Pc it) 

- ^ [ V2rq it) Pc it). (B14) 

^ Jx,t 

The modified propagators are given as 

(Gf )-i (r, t') = {-d^ + tm ir) - V^) d (t - r') 

+ (t - t') (BI5) 

on the imaginary axis and 

^-Mt,0 =((G«)-i '^^^^-rqit)Sit-t')l 

(B16) 

with modified retarded function 


(Gf) ^ it, t') = - (dt + Tc it) - V^) 6 it- t')+ 6 p it - t') 

(B17) 

on the real axis and 2x2 unit matrix 1. Thus, the bare 
masses have been replaced by in general time dependent 
terms that corresponds to become the usual self energy 
corrections of the large-TV theory. 

We express the vectors ^c,q,M terms of the com¬ 
ponent 4>c,q,M along the direction of the field = 


This is a crucial aspect of the theory demonstrating that 
the bare post-quench correlation function is affected by 
the pre-quench behavior due to memory effects of the 
bath. 

Next, we integrate over the transverse fields iTc^q and 
obtain 


S = 

+ 



^ [ (l>MiT)iGy) ^ iT,T')(l)MiT') 



(t>M (r) f] (it, t) y it) 


(B21) 


+ T / ^^9 

^ Jx,t 

+ *y / {ro,iPM (r) -f yPM (r) - tm (r) pm (r)) . 
+ 1 (iV - I) tr log gy + ^ (iV - 1) tr log(Gf )-i. 


The action is of order TV, allowing us to perform a saddle 
point analysis. Minimizing the action with respect to (j), 
p, and r yields the nine saddle point equations. Minimiz¬ 
ing with respect to the fields on the Matsubara segment 
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of the contour yields 

hM (t) = J {t') 

+ [ fj {iT, t) (j)q (t) 


PM (r) = (r) + 2 

TM (t) = fo,i + UiPM (t) 


iJtrlogtJ^^ (5trlog(Gr ) 


M\-l 


Stm (t) 


(B22) 


The minimization with respect to the quantum field on 
the Keldysh contour yields 


he (t) = -Tq (t) (j)q (t) + j dt'{Gf ) ^ (t, t') (j)c {t') 

+ j dt' (G~^)^ {t,t')(j)q {t')+ J (j)M (T)f]{iT,t) 


Tc {t) = roj + ^Pc (t) ■ (B23) 

Finally, from the minimization with respect to the clas¬ 
sical fields on the Keldysh contour follows that 


hq (t) 

Pq (t) 
rq (t) 


-Tq {t) (pe (t) + 

-^(pe (t) (pq (t) 

u . , 

-^Pq (t) ■ 


[ {G^)-^it,h)p>q it') 

Jt' 

i (Jtrlogt/F^ 

SVe (t) 

(B24) 


The natural saddle point solution for the quantum com¬ 
ponents respecting causality is rq = pq = (pq = 0 for 
hq = 0. We checked that this trivial solution is the only 
solution by an explicit analyses of the Heisenberg equa¬ 
tion of motion of the order parameter and the propaga¬ 
tors. Performing the functional derivatives and evaluat¬ 
ing them at vanishing quantum fields yields 

5trlogei-i 6g-^ 

6reit) V'^Sreit) 


Jtrlogg-^ 

Srq (t) 

(5trlog(Gf)-i 
Stm (t) 


= -Gf it,t)-G^ it,t)=0 

= (Gf)-i(r,r) 


Jtr log ^ 
SrM (t) 


(B25) 


We consider time-independent fields hjyj (t) = hi and 
he (t) = hf before and after the quench, respectively. 
Since the imaginary time evolution prior to the quench 
is in equilibrium we obtain time independent solutions 
for tm = Ti and (pM = <pi, yielding the usual equilibrium 
version of the large-equations. Eliminating pM yields 


In the last step we reintroduced spatial coordinates. The 
non-equilibrium large-A'^ equations of the classical com¬ 
ponents Ve, Pc, and (pe are nontrivial and time dependent. 
We obtain after the elimination of p: 

he= J dt'iGe)~^ it,t') (pc it')+(pM j fl{iT,t) 
r (t) = roj + (t) + (x, h x, t) . (B27) 

These are the self-consistent large-equations given in 
the main part of the paper. There we only used the 
more physically motivated notation hi = Am, (pi = (pM, 
and Ti = Tm to indicate that the variables on the Mat- 
subara branch refer to the initial field, order-parameter, 
and renormalized mass, respectively. Similarly we have 
hf = he, (pit) = (pc it) , and r (t) = rdt) for the corre¬ 
sponding variables after the quench. 


Appendix C: The fixed mass post-quench Green’s 
functions 

In this subsection we summarize the main steps in the 
derivation of the post-quench retarded and Keldysh func¬ 
tion of the system with constant masses r. The knowl¬ 
edge of this propagator is essential for any development 
of a perturbative approach to includes interactions. 

In order to determine the bare propagators, we start 
from the Heisenberg equation of motion of the field op¬ 
erator after the quench: 

nOO 

{dt + roj + q^) ^ iq, t)= dsr] it - s) ip (g, s) 

Jo 

+ aiq,t) + hiq,t) (Cl) 

with external field h iq, t) that is coupled to the order 
parameter. The source operator Siq,t) is given by 

“( 9.0 = -/ dsvit-s)ipiq,s). (C2) 

J —00 

It is useful to express H(t) in therms of the initial bath- 
operators Xj = Xj iq, t = 0 ) and P° = Pj iq, t = 0 ): 

“ i) = - XI 

j \ J / 

We solve the Heisenberg equation of motion via Laplace 
transform 

nOO 

ipiq,u;)= dte<-+'^^>ipiq,t) (C3) 

Jo 

with boundary condition tp iq,t = 0) = ipg iq) and 

dt<p iq, t = 0) = ttq iq). It follows 

ipiq,uj) = Fiq,oj)gfiq,uj) (C4) 


hhi = rM(pM 

Tm = ro,i + '^(pii+UiGd ix,T-,x,T). (B26) 


with force-operator 

F iq, w) = TTo iq) - iujipQ iq) -f H (q, w) -f h (g, w) (C5) 
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and bare retarded post-quench Green’s function 


gf (9,w) 


1 

- tqj - q'^ + r]{u})' 


(C6) 


system should of course have the same formal structure 
as Eq. (C9), i.e.: 

gf {g, w, w') = M {q, w, w') gf {q, w) gf {q, uj') (Cll) 


To check that this is indeed the correct result for the 
bare retarded Green’s function one expresses the back 
transform ip {q, t) of p {q, oj) in terms of gf {k, w): 

cp {q, t) = (tto (g) + ipQ (q) dt) gf {q, t) . (G7) 

pOO 

+ (3(g,s)+h(g,s))gf (fc,t-s)ds 

Jo 

Inserting this result into the definition Eq. (37) of the re¬ 
tarded function does indeed yield Eq. (C6) for zero exter¬ 
nal field h (g, t). Thus, for the bare retarded post-quench 
Green’s function follows the same result as in equilib¬ 
rium, which only depends on the difference between the 
two time arguments. As we will see below, this behavior 
of the retarded function will not carry over to case where 
we include interactions. 

Next we consider the bare post-quench Keldysh func¬ 
tion gf (g, t, t'). This function will depend on both time 
scales already for a non-interacting system. We deter¬ 
mine the Gf (g, t, t') using the same approach as for the 
retarded function, i.e. insert the solution Eq. (C7) into 
the dehnition Eq. (38) of the Keldysh function. Here, it 
is convenient to consider the double Laplace transform: 

gf (g, a;, cc') = r dtdt'gf (g, t, t') e<-+*+“V*'), (C8) 
Jo 

with = to iO'^. From our solution Eq. (C4) follows 
gf (g, w, w') = M (g, w, w) gf (g, w) gf (g, w'), (G9) 

with memory function 


The subscript i indicates that we are considering a system 
in equilibrium that is governed by the Hamiltonian Hi 
with bare mass ro,i. The key insight is that the memory 
function M(g,w,w') must be the same function as in 
Eq. (G9) as we have to determine the expectation value 
of the same operators F (g, oj) with respect to the same 
state. 

In equilibrium, the Keldysh function only depends on 
the difference of the time arguments. For the double 
Laplace transform this implies 


gf (g,w,a;') 


, gf jg,^) + gf 

UJ + uj' + 10 + 


(G12) 


where gf (g, uj) is the Laplace transform of the bare 
Keldysh function of the initial state in equilibrium: 

pOO 

gf{q:‘^)= (g,t). (G13) 

Jo 

In equilibrium we can then use the fluctuation-dissipation 
relation 


gf (q, f = ( 2 ^^) 


to determine this function. Thus, we obtain for the mem¬ 
ory function 


M (g, w, uj') = i 


gf{q,uj)+gf{q,uj') 
UJ -\-uj' + i0+ 

-1 I 


X gf (q, w) gf (g, uj') 


(CIS) 


M{q,uj,uj') = -i{[Fa{q,uj) ,Fa{-q,uj')]_^_) , (CIO) 

given by the force-force correlation function. M (g, w, uj') 
can be obtained from the expectation values of F (g, uj) 
as given in Eq. (C5). To proceed we note that the oper¬ 
ators that enter the memory function are Pq (g), (g) 

etc., i.e. they are for system and bath variables at t = 0. 
Thus, all relevant expectation values can be evaluated 
in equilibrium prior to the quench. The direct evalua¬ 
tion of these expectation values is straightforward but 
somewhat tedious. In particular one finds that numer¬ 
ous expectation values diverge in the limit of an infinite 
cut-off uJc, divergences that cancel if one combines all 
terms that contribute to the memory function. To avoid 
these complications we present a signihcantly easier ap¬ 
proach to this problem. We explicitly checked that both 
approaches lead to the same result. We stress again that 
all expectation values that enter the memory function 
can equally be determined prior to the quench when the 
system is still in equilibrium. We therefore consider a 
system without quench and with initial Hamiltonian Hi 
for all times. The Keldysh function of this equilibrium 


As mentioned, we obtained the same result using 
a straightforward evaluation of the definition of the 
Keldysh function, inserting the solution Eq. (C7), and 
evaluating all expectation values explicitly. As limiting 
cases, Eq. (C15) includes the solution of Ref. 12 for a clas¬ 
sical system with ohmic bath and Ref. 81 for a classical 
system with colored noise. In addition it reproduces the 
findings of Ref. 86 for a non-interacting quantum system 
without bath. 

As discussed in the context of the scaling behavior, 
the limit where we consider large distances Srgi from the 
critical point prior to the quench is of particular interest. 
In this deep quench limit follows: 


M{q,uj,uj') 


,coth(^)(5?7(a;) -I- coth{^)6r]{uj') 
UJ + uj' + i0+ 


(G16) 


For times t,t' > U = the memory function 

is always governed by this deep-quench limit, even if dvi 
is nominally not large. This enables us to introduce a 
scaling form for and to expect universal behavior af¬ 
ter the quench. The justification to use the deep-quench 
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limit in this paper is however not, that we concentrate 
one time-scales larger than ti, but rather that the even 
for a small quench the spatial correlation length collapses 
and the bar Green’s function become completely local in 
space, corresponding to a large effective mass at time 
t = , independent of the quench protocol. 


Appendix D: Long time limit of the Green’s 
functions 

We can extend the analyses for the equal time Keldysh 
function to different time arguments t,t', but still in the 
limit 

^ ^ < 1. (Dl) 

We use the long time expansion for (p ~ peq — 
It dsG^^it - s)r{s)<peqis) to obtain: 


The same procedure leads for the second integral in 
Eq. (D2) : 

/■*' 

- J dsG^^it' - s)r(s)G^(s - t) 

^ -rit') I (D6) 

For r{t) ^ r{t') « r(G) with G = {t + t')/2 this yields to 
the result given in Eq. (103) in the main text. In analogy 
to the long time limit of G^{q, t, t) we can derive a simi¬ 
lar behavior for the retarded Greens function G^{q, t, t'). 
It follows with the long time expansion for cp and the 
definition of G^ in Eq. (37) that 

G^{t, t') = -ie{t - 

+ ie{t-t') dsG^^it- s)r{s){[(p{s),ip{t')]-) 

+ ds{[ip{t), vj(s)]_)r(s)Gfq(t' - s)^ 


G^{q,t,t') = - i{[p{t),p{t')]+) 


=Geq{q,t - t') - [ dsG'fq(t-s)r(s)G^(s,t') 
Jt^ 

- f dsG^^it' - s)r(s)Gfq(s - t) (D2) 


To evaluate the s-integration we express G^ by its 
Laplace-transformation and G^ by its Fourier transfor¬ 
mation: 


-J rfsGfq(t- s)r(s)G^(s,t') 


= — I 


"^^ImGfq(w)G^(w') f dsr(s)e-“(‘-")-*“'(*-*') 


2 Tr^ 


(D3) 


=Gfq(t - t') - £ dsG^^it - s)ris)G^^{s - t'). (D7) 

We express the equilibrium G^ by their Laplace- 
transformation and obtain 

(5Gf (t, t') =- dsG^qit - s)r(s)Gfq(s - t') 

= j ^^^^ImG'fq(w)ImGfq(w') 

r ^ 

X / (D8) 

Jt' 

uj and uj' are of order of the typical mode-frequency 
q“^7^/2. In the limit where the time difference t — t' 
is small, but t, t' are both large compared to the mode- 
frequency the s-integral can be approximated by 


Here, w and w' are of order of the typical mode-frequency 
In the limit where the time difference t — t' is 
small, but t, t' are both large compared to the mode- 
frequency the s-integral can be approximated by 


dsr(s) 




+ 1 ( 
\ 2 ) i{uj — uj') V 




')) (D9) 


/■* p-ibj' (t-t') 

I dsr(s)e-“(*-")-“ ~ r (t) ^ (D4) 


In a next step, we use Kramers-Kronig relation to pre¬ 
form one frequency-integration. This finally yields 

-|‘dsGfq(t-s)r(s)Gfq(s,t') 

= - r{t) I ^ReG3« (D5) 


This yields to the result given in Eq. (104) 
SG^{t t') ~ -ir ^ 


(' 




= —4dr 


2*r I ^ReGfq(a2')ImGfq(a;')e-“'(‘-‘') 

* + G^{t-t'). (DIO) 
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